Can Induced 6 Vacua be Created in Heavy Ion Collisions? 
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The development of the early Universe is a remarkable laboratory for the study of most nontrivial 
properties of particle physics. What is more remarkable is the fact that these phenomena at the 
QCD scale can be, in principle, experimentally tested in heavy ion collisions. We expect that, in 
general, an arbitrary induced 6 vacuum state d^'"'')) would be created in heavy ion collisions, similar 
to the creation of the disoriented chiral condensate with an arbitrary isospin direction. It should be 
a large domain with a wrong 9^"'^ 7^ orientation which will mimic the physics of the early universe 
immediately following the QCD phase transition when it is believed that the fundamental parameter, 
Qfund _^ Q -y^g ^gg^ ^j^jg jjg^ numerically in a simple model where we study the evolution of the 
phases of the chiral condensates in QCD with two quark flavors with non-zero ^'""^-parameter. We 
see the formation of a non-zero ^^"''^-vacuum with the formation time of the order of 10"^'^ seconds. 



1. The realization of the Relativistic Heavy Ion Col- 
lider (RHIC) at Brookhaven opens up some exciting 
doors in fundamental physics research. The possibility 
of producing and studying quark-gluon plasma is the 
most well publicized research area. But it is certainly not 
the only new field that RHIC makes possible. This let- 
ter is concerned with the possible creation of non-trivial 
^""^-vacua inside the high temperature fireball created at 
RHIC. 

The idea is very similar to the creation of the Disori- 
ented Chiral Condensate (DCC) in heavy ion coUisons 
B, P], (see also a nice review 0] for a discussion 
of DCC as an example of an out of equilibrium phase 
transition). DCC refers to regions of space (interior) in 
which the chiral condensate points in a different direction 
from that of the ground state (exterior), and separated 
from the latter by a hot shell of debris. As we shall see 
in a moment, for both cases (DCC and |6'""')-state) the 
difference in energy between a created state and the low- 
est energy state is proportional to the small parameter 
TOq and negligible at high temperature. Therefore, ener- 
getically an arbitrary |0*"'') can be formed which is the 
crucial point for what follows. 

Before going into details, we would like to recall some 
general properties of the DCC, which (hopefully) can be 
produced and seen at RHIC, with an emphasis on the 
analogy between the DCC and a misaligned |6')-state. 
If the cooling process is very rapid and, therefore, the 
system is out of equilibrium, there will be a large size 
of the correlated region in which the vacuum condensate 
orientation mismatches its zero temperature value. The 
absolute value of the chiral condensate right after the 
phase transition is expected to be close to its final (zero 
temperature) magnitude. However the vacuum direction 
of the formed condensate is still misaligned since it takes 
a longer time for the vacuum orientation to relax due to 
the small free energy difference '^ rriq between the formed 
and true vacuum states. 

Due to the fact that the absolute value of the chiral 
condensate right after the phase transition is expected 



to be close to the zero temperature magnitude we can 
parametrize Goldstone fields by matrix U in the following 
way: 

U = e'^("-^) , C/C/+ - 1 , {^\^'j,) = -\{^L^u)\ U,, 

The energy density of the DCC is determined by the mass 
term: 

E^ = — rr(Mt/-|-Mt[/t) = -2m|(*1')|cos(</>) (1) 

where we put 7n„ = tti^ = vn for simplicity. Eq.(|l|) im- 
plies that any (f) ^ 0{mod 27r) is not a stable vacuum 
state because -q^\^^o ^ 0, i.e. the vacuum is mis- 
aligned. On the other hand, the energy difference be- 
tween the misaligned state and true vacuum (0 = 0) is 
small and proportional to rriq. Therefore, the probability 
to create a state with an arbitrary at high temperature 
T ^ Tc is proportional to exp[~V{E^ — Eo)/T], where 
V is 3D volume, and depends on cf) only very weakly, i.e. 
(/) is a quasi-flat direction. Right after the phase transi- 
tion when (^^) becomes non-zero, the pion field begins 
to roll toward (p = 0, and of course overshoots = 0. 
Thereafter, oscillates. One should expect the coherent 
oscillations of the tt meson field which would correspond 
to a zero-momentum condensate of pions. This is ex- 
actly what was found in g]. Eventually these classical 
oscillations produce real tt mesons which hopefully can 
be observed. 

Now we turn to our main point when the U{1)a phase 
of the disoriented chiral condensate is also non-zero and, 
therefore, the |6')-vacuum state could be formedn. The 



*From now on we omit the label "ind" for the induced 6. 
We hope S'"'* will not be confused with g-f""'' which is zero in 
our world and which can not be changed. The simplest way 
to visualize S'"'* is to assume that right after the QCD phase 
transition the flavor singlet phase of the chiral condensate 



production of non-trivial 0-vacua would occur in much 
the same way as discussed above. The new element is 
that in addition to chiral fields differing from their true 
vacuum values the ^-parameter of QCD, which is zero in 
the real world, becomes effectively nonvanishing in the 
macroscopically large region. 

In this letter we show that in a simplified numerical 
model non-trivial 0-vacua can be realized. While this is 
a simplified model which differs somewhat from the real 
physics associated with high energy ion collisions, there is 
no reason to believe that the non-trivial 0-vacua cannot 
be realized at RHIC. Our results also give a very rough 
estimate of the time it takes for these non-trivial 0- vacua 
to be formed. In order to have a hope of observing them 
they must form within the time that the central region 
of the fireball is isolated from the true vacuum. 

2. To take into account the U{1)a phase associated 
with |0)-vacua we choose the matrix Uij in the form U = 

diag (e*"^') with J2i=i 4'i in general non-zero. The energy 
density of the misaligned vacuum is determined in this 
case by the following low-energy potential 0,0]: 
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Z = 0,l,...,iVc-l, 



(2) 



where E = {baJ{Z2'K)G'^) - IQ^'^GeV^ is much larger 
than M, = -m^(^*} - IQ-^GeV^. Here h = ll/iN^ - 
2/3Nf, rUg - 5MeV, (^^) ~ -{240MeV)^. The crucial 
point is that the 9 parameter appears only in the com- 
bination, '^ (/)i — 9. This is a direct consequence of the 
transformation properties of the chiral fields under U{1)a 
rotations. To convince the reader that (0) does indeed 
represent the anomalous effective low energy Lagrangian, 
three of its most salient features are listed below: 

i) Eq. (0) correctly reproduces the Witten-Di Vecchia- 
Veneziano effective chiral Lagrangian jg] in the large Nc 
limit; 

ii) it reproduces the anomalous conformal and chiral 
Ward identities of QCD; 

iii) it reproduces the known dependence in 9 {ie. 2tt 
periodicity of observables) [g[. As mentioned above, in 
the large Nc limit the effective Lagrangian (g) takes the 
Witten-Di Vecchia-Veneziano quadratic form. We would 
like to remark here that the exact form is not impor- 
tant for our present purposes as long as the 0-parameter 



appears in combination '^ cjji — 9 and the parameter 
E > M,. 

The most important difference between Eqs. (0) and 
(0) is the presence of the parametrically large term ^ 
E :^ m, 1(^^)1 in the expression for energy (^), describ- 
ing the U{1)a phase of the disoriented chiral condensate. 
This term does not go away in the chiral limit and pro- 
vides a non-zero mass for the t]' meson which is expressed 
in terms of the parameter, E. It was for exactly this rea- 
son that it was thought until recently H that the non- 
trivial 0- vacua would involve too large an energy cost to 
be produced because of the large parameter associated 
with the 9 parameter. 

The key point is the following. For arbitrary phases 4>i 
the energy of a misaligned state differs by a huge amount 
^ E from the vacuum energy. Therefore, apparently 
there are no quasi- flat (~ niq) directions along (jji coor- 
dinates, which would lead to the long wavelength oscil- 
lations with production of a large size domain. However, 
when the relevant combination {^.^4>i — 9) from Eq.(||) 
is close by an amount '^ O^niq) to its vacuum value, a 
Boltzmann suppression due to the term ~ ii^ is absent, 
and an arbitrary misaligned |0)-state can be formed. 

Indeed, in the limit Mi <C E, because of the large pa- 
rameter in the first term the vacuum state that is favored 
for non-zero 9 is the solution (pi « 9/Nf. Substituting 
this back in to the potential we reproduce the well-known 
result [pj for the vacuum energy as a function of 9. 



V{9) K-E-J2 Mi cos{9/Nf) 



(3) 



This shows that the energy cost of creating a non-trivial 
0-vacuum goes like the much smaller parameter Mi in 
agreement with a general theorem that the 9 dependence 
appears only in combination with niq and goes away in 
the chiral limit. 

At this point we can apply the same philosophy as for 
DCC. The chiral fieldsH (pi are allowed to take random 
values and after the phase transition begin to roll toward 
the true solution (j>i « 9/Nf and of course overshoot it. 
The situation is very similar to what was described for the 
DCC with the only difference that in general we expect 
an arbitrary |0)-disoriented state to be created in heavy 
ion collisions, not necessarily the \9 = 0) state. The 
difference in energy between these states is proportional 
to niq from (0). 

To be more quantitative, one should study the evo- 
lution of the chiral fields with time. For this paper we 
assume the following situation. The rapid expansion of 



is non-zero in a macroscopically large domain. This phase 
is identified with 0'^"'', This identification is a direct conse- 
quence of the transformation properties of the fundamental 
QCD lagrangian under U{1)a rotations by which the chiral 
singlet phase can be rotated away at the cost of introducing 



^If 6* 7^ 0, the Goldstone fields are not exactly the pseu- 
doscalar fields, but rather are mixed with the scalars; the 
mixing angle between the singlet and octet combinations also 
depends on 6, see |7| for detail. 



the high energy shell leaves behind an effectively zero 
temperature region in the interior which is isolated from 
the true vacuum. The high temperature non-equilibrium 
evolution is very suddenly stopped, or "quenched" , leav- 
ing the interior region in a non-equilibrium initial state 
that then begins to evolve according to (almost) zero 
temperature Lagrangian dynamics. Starting from an ini- 
tial non-equilibrium state we can study the behavior of 
the chiral fields using the zero temperature equations of 
motion. The equations of motions are non-linear and 
cannot be solved analytically but we can solve them nu- 
merically in order to determine the behavior of the fields. 
The production of a non-trivial 0- vacuum is indicated by 
the fact that the chiral fields relax to constant and equal 
non-zero values on a time scale over which spatial oscil- 
lations of the fields vanish. This fact confirms the self 
consistency of our approach. First we assumed that the 
singlet phase of the chiral condensate is non-zero on a 
macroscopically large domain. We then calculated the 
evolution of the chiral phases in this background. Finally 
we closed the circle by realizing that the chiral fields relax 
to constant, equal, non-zero values as mentioned above. 
The formation of a non-perturbative condensate is also 
supported by observation of the phenomenon of coarsen- 
ing (see below) and by a test of volume-independence of 
our results. 

3. The equations of motion for the phases of the chiral 
condensate with two quark fiavors consists of two coupled 
second order nonlinear partial differential equations: 



^,-v20, + 70, + -^y(0,,0) 



i = l,2 (4) 



where V^ is a three dimensional spatial derivative and 
the potential is given in (0). Emission of pions and ex- 
pansion of the domain will contribute to the damping, 7, 
as might other processes. We do not know exactly how 
they would contribute but we simulate these unknown 
effects by including a damping term with a reasonable 
value for the damping constant, 7 ~ Aqcd ~ 200 MeV. 

The initial data for each of the chiral fields 0i is cho- 
sen on a 3D grid of 16'^ points. The initial data con- 
sisted of random values of 4>i and (pi — 0. The initial 
data was evolved in time steps using a Two-Step Adams- 
Bashforth-Moulton Predictor-Corrector method for each 
grid point with the spatial laplacian approximated at 
each grid point using a finite difference method. We used 
periodic spatial boundary conditions. 

The grid spacing was determined by the length of a 
side of the spatial grid which was varied in order to vary 
the volume. The size of the time step between succes- 
sive spatial grids was much smaller than the spatial grid 
spacing and was fixed at about lO~^MeV~^. 

We evolved the data for 8000 time steps and then ap- 
plied a Fast Fourier Transform [g| to the spatial data at 
evenly spaced time steps. We then binned the data in 
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Time(1.6x10-tMeV-i) 
FIG. 1. |(/f>fc| is shown for various |fc| as a function of time. 
Notice that the zero mode settles down to (pi ~ d/Nf. 



small increments of the magnitude of the wave vector in 
order to obtain the angular averaged power spectrum. 

This procedure was carried about for different volumes. 
In all cases the results were qualitatively the same. We 
saw an initial growth of long wavelength modes as in 
B and subsequent damped oscillation of all modes. The 
|fc| = modes oscillate and approach the equilibrium val- 
ues of the fields. They exhibit this behaviour in the same 
time frame in which the Fourier coefficients of the modes 
with non-zero wave vectors fall to a tiny fraction of the 
zero mode coefficient. This qualitative behavior occurs 
for different total volumes and grid sizes suggesting that 
this behavior is not due to finite size effects. 

We should note that our |/j| = mode is really only 
a quasizero mode as it is obtained in a finite spatial vol- 
ume with periodic boundary conditions. However, our 
quasizero mode approaches the same value irrespective 
of the total spatial volume indicating that this really is a 
condensate. If it were not we would expect the value of 
the coefficient to decrease when the volume of the system 
increases. 

The evolution of the Fourier modes of the (j> fields is 
shown in Fig.H for the specific case oi 6 = 27r/16 and a 
spatial grid of 10 fm on a side. The initial values of the (p 
fields were randomly chosen within the range —77r/16 to 
Ttt/IG. The zero mode clearly settles down to a non-zero 
value. All higher momentum modes vanish extremely 
rapidly and are negligible long before the zero mode set- 
tles down to its equilibrium value. 

The instantaneous distribution of Fourier modes for 
the evolution above is shown in FigJ^ at a few different 
times. This graph clearly shows the amplification of the 
zero mode as time increases. This phenomenon of coars- 
ening and the formation of a nonperturbative condensate 
is very similar to earlier discussions in ref. y. 

In Fig|3| we plot |0fc| as a function of time for three 
different volumes. We chose 9 = 7r/16 and the volumes 
{8fm)^, {16fm)^, and (32/m)'^. For each volume, we 
plotted the zero mode and the same non-zero mode. No- 
tice that the zero mode is independent of the volume of 
the system, while the magnitude of the non-zero mode 
decreases with increasing volume. This is the signature 
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Magnitude of k (arbitrary units) 
FIG. 2. Above we demonstrate that the system exhibits 
the coarsening phenomenon. The data was sampled at three 
times within the first 1000 timesteps of the evolution. 




Time (1 .6x1 0-tMeV-i) 
FIG. 3. The zero mode and a non-zero mode are shown 
as a function of time for three different volumes. The heav- 
ily dashed line represents the smallest volume, the medium 
dashed line the middle volume, and the solid line represents 
the largest volume. This graph illustrates the volume inde- 
pendence of the zero mode. 



correlations suggested in ref. [Q. However, one should 
expect that the signal will be considerably (if not com- 
pletely) washed out by the re-scattering of the pions and 
their interactions in the final states, which mimic true 
CP-odd effects. In practice it is quite difficult to over- 
come the problem of separating a true CP violation from 
its simulation due to the final state interactions. A more 
promising direction is to to look at r]{r]') -^ tttt [p^ 11 
decays which are strongly forbidden in our world, but 
nevertheless will be of order one ii 6 ^ 0. Indeed, one 
can show ffll that 



T{r, 



2 1- 



4w4(fl) 



\{0\qq\0) 



0.5 MeV{sm-f, 



f3 



(5) 



and therefore the effect could be quite noticeable. We 
should remark here that all masses in the 6 vacuum state 
are shifted in comparison with their values at 6* = M. 
Formula (pi) reduces to the corresponding expression of 
reference |]l2| in the limit 6 ^> 0. An analogous calcula- 
tion for rj' — > TTTT leads to a similar numerical estimation 
for r{r]' -^ tttt) ^ 2 MeV{sm |)^ which is almost an 
order of magnitude larger than the full width of the 77' 
meson in the 9 — world. More importantly, this width 
is exclusively due to the CP odd decay r(?7' -^ tttt). 



that a real condensate has been formed. 

For a total volume of (10/m)"^ and 9 — 27r/16 the time 
for relaxation from the initial non-equilibrium state fol- 
lowing the quench to the non-trivial 0- vacuum is approx- 
imately 0.064 MeV^i « 4 x lO^^^s. This seems to be of 
the same order of magnitude as the time which we might 
expect the central region of the fireball to be isolated 
from the usual vacuum. The volume we have used is just 
at the upper limit of what we would expect at RHIC. As 
well we have not attempted to correctly account for the 
collision geometry at RHIC either. However, our sim- 
plified calculation suggests the possibility of producing 
non-trivial 0-vacua. 

4. We have shown through a numerical calculation of 
the zero temperature equations of motion with a non- 
zero induced 9 parameter that the chiral fields (f> after a 
quench from high temperature go to a spatially constant 
non-zero value related to the 6'-parameter. This occurs 
on a time scale of the order of 10"^'^ seconds. The fact 
that all other non-zero modes fall to negligible values long 
before indicates that we have formed a condensate or a 
non-trivial \9) vacuum state. 

The most intriguing question is: "What is the signa- 
ture of the produced 9 state?" First of all, due to the 
fact that the 0-vacua are odd under charge conjugation 
times parity, CP, their decay must produce some CP odd 
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